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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM NO. 1195 


ON THE SOUND FEED) OF A ROTATING RR0RELI2R* 
By L. Gutin. 


The sound field of a rotating propeller Is treated theoretically 
on the Basis of aerodynamic .principles. For the lower harmonics, 
the directional characteristics and the radiated sound energy are 
determined and are in conformity with existing experimental results. 


1. INTRODUCTION 

A rotating propeller produces per iodic “disturbances of the 
medium vhioh cause a sound of low frequency. The fundamental tone of 
this Bound equals the product of the number of blades and the numb er 
of revolutions per unit time, since each .spatial -configuration of 
the propeller is repeated with exactly the same frequency. . The ■ 
present work a-ime at a theoretical investigation of the sound field 
produced by the propeller. First's short report on earlier investiga- 
tions in this field will be given. 

The directional properties of propeller sound were first investi- 
gated theoretically by Iynam and Webb,-*-' They suggested two hypotheses 
on the acoustic action of the propeller. Both hypotheses were based on 
actually arbitrary assumptions and led to directional characteristics 
(cf. figs; 1 and 2) “which are not in agreement with the observed ones 
(cf. figs. 3 and 4). 

Hart® attempted to develop a theory which was supposed to be free 
of arbitrary assumptions; however, he tacitly made the same assumption 
as Iynam and Webb in their second hypothesis: namely, that each 


*"Uber das Schallf eld einer rotierenden Iaiftschraube . ” ' Physikalisehe 
Zeitschrift der Sowjetunion, Band 9, Heft 1, 193^, PP* 57-71* 

^Tynaa, E, ' J, and Webb, H. A.t The Emission of Sound by Airscrews, 

R, & M. No. 624, 1919. (The theories of Iynam and Webb are known to us 
only through the treatment of Baris and Kemp since the original report 
was unfortunately not accessible.) 

-SHart, M. D. : The Aeroplane as a Source of Sound, R. & M. No. 1310, 
1929. 
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disturbed element in the propeller plane acts as a simple nondiroctional 
source; accordingly , he obtained ’the ~saae' ‘result. 


■ ■ : ■ r;.\ .o- :• ,• • ‘ a 4, 

Experimental investigations were undertaken by Paris and Kemp . . 
Both experimented with a two -blade propeller; its diameter was 
4.5 meters, the number of revolutions 13*9 per second. It is true 
that a notioeable discrepancy exists between their results in quantitative 
respect; • essentially., .however ,. they . agree . Both directional character- 
istics show an . asymmetry with, respect to the rotational' plane, a 
principal maximum at ; about . $ =i 115°>' and a secondary maximum at 
d a 40 ». This peculiar a eymmejt^.- .contradicts both hypotheses by 
Lynam and Webb.. 


Paris obtained, by a combination of the two hypotheses, a 
directional characteristic (of . fig, 6) which reproduces the main 
outlines of the experimental- curves; however, his" method lacks 
physical foundation, and. the number of arbitrary assumptions is too 
large. . \ . . .. V . 

The hypotheses mentioned above did not yield any information on 
the sound output of the propeller* This problem was experimentally 
examined by Kemp. He obtained the values 17.8 watts for the first 
harmonic and 7*4 watts for the second'. ® . 

, ’ . * t 

In the following, a theory based on aerodynamic principles shall 
be developed, which determines quantitatively the sound field of a • 
propeller. - ■ • 


2. SOMES AERODYNAMIC 'RESULTS 


A few facts concerning the. action of the propeller,, which are 
known from aerodynamics, shall" be briefly quoted. 


The cross section of a propeller blade looks very similar to 
that of a- wing (of; fig. 7) j. the action 'of the propeller is based on ' ‘ 
this fact* Motion of a wing element relative, to the medium c.au-es, 
exactly as for a wing, a pressure increase oh the concave side and 


Paris , B. T. Phil*,; vpl ; ‘ 13 ,' : no s'Jjjp ‘ 1932 . • , <: ‘ 

ietitp, 'Si' F./Proc ."Physical SoSH (‘Bdhdpnl^^ql. ',44/.. pi". 2, .. _p.* 151'. 


In .the report by Kemp values half as large as these ore given 
by mistake. Compare ‘Paris; 'E. T.y'ph'ili' M^Pvoi-.' 15; no i 60, 1933* 


6 Compare with this paragraph Mises, R. 7.5 The Theory of Plight, 

pp. 116-121* 
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a pressure decrease: thsai; cdnvex side'. -The: resultant ’‘iir .f or ce may 
Id© ■ separated intO.-two components.' 'of >rM<Sh oils has the direction of/ 
the propeller, axis and indeed is a. thrust in the direction, of '.the 
forward motion, the other . one is opposed: to the rotation and represents 
a drag force, .r The thrust, forces - of-- single elements add up to the- 
total thrust/f orce of the.propelle'r, the drag forces result 'in a 
torque opposed to the, torque of th& motor. ’ '. ' ! 

A : few-i%03?tant quantitative results which will. be., used subsequently 
are derived from the so-called momentum theory of the propeller. ; .The 
propeller is assumed to "be at rest and the air moving towards it from 
the front, opposed to the flight direction, .vith the velocity- V. The 
effect of the ‘pr.<^iler/^n consists 'in' an : acceleratioh- of the,.', . 
approaching mass of , air which' finally assumes at a certain distance / 
behind the -propeller the .velocity' V + w. Let p be .the density of ' 
the medium, 'P /the magnitude of the; thrust force, 1 ’ S .'the' area /, , * . 
described by the propeller, * W ' the power supplied to the propeller, . Jftj 
the part of the power converted into translational energy. Then the 
following relations hold: . .. . - 



By eliminating v .from both equations one obtains an expression 
for the dependence of- the thrust on the "power. For static conditions 
(V s* 0) there results: • / ■ .... 


The obvious relation ■, 








shall be mentioned. ' :(M = torque of .the .motor eves angular velocity' . • 
of the rotation.) / ■' * 
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3. D3OTRBANC3B FORCES IN THE PBOFEIXER PLANE 

Two foroesact- on each element- of the .propeller : a thrust and 

a drag force* From the well-known theorem of mechanics, It follows 
that each element -exerts forces of equal magnitude and in the 
opposite direction upon the medium* ' The points of application of these 
forces may he imagined concentrated in one plane since the axial 
dimension? of the propeller, are very small in comparison with the 
wave lengths of the. first harmonics j this plane is subsequently denoted 
as "the plane of rotation*. 

A propeller element will he considered# the distance of which from 
the axis equals R$ let dR he its radial length, a ' its width, 
measured in the projection on the rotational plane . Let the forces 
exerted upon.. the medium hy the element he A(R) dR (in the- direction 
of the axis, opposed to the flight direction) and B(R) dR (in the 
direction of the rotation). 

It is clear that 

fl E 0 - { 1*0 

n A{R) OR m <3P * ? 

I/O (lo 


1*0 

to 


n representing the number of blades, R(j the length of the blade. 

Let it first he assumed that the forces are uniformly distributed 
over the width of the element* 

In the element R cLR 09 of the rotational plane, forces 
B cL9 E <3*0 

A(R) dR— — and B(R) dR— act on the medium during the time 

£L 8, 

interval in which this element is covered hy the projection of the 
propeller element. If -the overlapping starts at the time t «= O' it 
will he ended at the time t « t ** to return after ' t * T = g*. 

These periodical forcos may he developed in Fourier series 
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FlCt) 


A(R}~~dR (0 < t < t) 
& 


* / Am cos (irnat - 6 m ) + A Q 

o . . c t < t <t) fe 


CO <t < t) ^-SL : - 
Fo(t) = % «* / Bju, oos (renat - t^) + Bq 

c o . . (t <t <t). . 


On© obtains 


Am = ~A(E)| sin (maj) cffi 00 
B m = sin SR d£ 



(3) 


In a second area element R dR dS, shifted with respect to 
the first by the angle 6 in the rotational direction, there act 
periodical forces of the sane magnitude but retarded by the time * 
The corresponding Fourier developments are 




Aja cos (mnat - mn9 — e m) 


+ A o 


w 

> Bm cos (mnat - mn0 - %) + Bq 


m=l 


For the first harmonics 2121 - i B usually small (the blade parts 

lying close to the center where R is small are eliminated since 
they make almost no contribution to the air forces), and one may equate 


sin 




ss 
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One then obtains 


• Bm 


2 A(R) dB 40 
§ B(R) as 40 



W 


It shall be shown now that these expressions are vali4 also when 
the air forces are not uniformly 4istributo4 over the blade width. 
Using for instance for the distribution of thrust the relation 


A(R) as 


l a 

A(R) as f(a) ~ 

do 


with 


f(s) ds * a 


(s is counted from the leading edge of the blade), one obtains for 
the Fourier coefficients the expression 


Ajn o §A(R)§ aR d0 


* T f(^f)a w i2DKt/D 4t 

Jo w 


For the first harmonic 
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and hence 


Am * §A(B) d B cL0 

Jt 


For e m one obtains 


6 m = arc 



sin 2nw~ 

T r 

<arc tan ■ » — ■■■ ■■ = 2nw— 

oos 2mrt- * 

a? 


A more accurate estimation gives 


ll 0 


( ' 2iturr 

for — 

f at cos 2 2mrt-^ + sin 2 2mrt^ 



with and tg, representing certain mean values. One can name 
distributions for which t^_ = 0 and tg « t, but for the distributicns 
occurring in practice the difference between t£ and tg is much 
smaller. 
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4. THE SOUND FIELD 


Tbs point of origin of coordinates is assumed at the propeller 
center, the y-axis and z -axis in the rotational plane (cf. fig. 8)* 

The x-axis is assumed to coincide with the flight direction. In the 
rotational plane polar coordinates are introduced, with tho polar 
angle being counted from the y-axie, in the direction of the rotation. 
Let it first be assumed that the rotation, observed from the flight 
direction, is counter clockwise. The forcos acting on the medium in 
an element E dE d0 are (for the first harmonics) 


X = dR d0 

Y = -S-B(R) sin 0 eiftost-mne-na) ^ d0 

Z = S b(r) cos 6 e i ( lcct " lnn0 *'%) as 




J 


(5) 


If one assumes, for reasons of simplification, that at the time 
t as 0 the center line of one of the blades coincides with the positive 
y-axis, 6m and t&h will be at any rate smaller than > therefore 

small, at least for the lower harmonics. 

The velocity potential oroducod by a concentrated force with 
the components X, Y, Z is ° 



( 6 ) 


8 Compare Lamb, H. : Textbook of Hydrodynamics, Teubner 1931, P* 567* 
A concentrated forco Fewest ie t therefore, equivalent to an acoustic 
double source of the strength 


kep* 
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If one inserts the expressions for the force components and 
notes that 



( % V are the angles that the radius vector r makes with the axis) , 
one obtains for the total velocity potential 


9 


in 


f) E 0 


E n ft 2rt 

[A(E)e ltkcVlra ' 3 - s i»> cob i 

Jo L 


lw s pok Jo 

+ B(E)6 1 ( kot - ,me - , >»>)(oos X sin 9 - cos v cos 0)J i~ aE as ( 7 ) 


The point at which the force acts is assumed to lie in the 
xy ■» plane which can always "be attained hy a suitable choioe of the 
system of coordinates j ,z in the rotational plane; let it further 
be assumed that r is large . Then it follows that : 


cos X = sin dj cos V a 0 


r = rp - E cos. @ sin d 
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e iKR sin 0 oos 6 


After substituting into equation (7) &nd carrying out the integration 
with respeot to 9 one obtains 


9 


jttPn e ll£ ( ot ~ r l) j 
2irpc r^ I 


E 0 


<3R 


A(R)e“ ie m cos d sin $) 


-in„ Jmn~l (KR sin d) + (l« sin d) 

2 


- B(R)e" iTte sin d 


i^^n e lk(ot - ri? 
2rtpc r d 


SOr 


i/O 


-A(R)e"' iem cos d 


+ ~B(R)e" i,,m J^hR sin d) dR 


(8) 




e iz cos (p-imcp dq) = 2 rti m j m ( z ) i 


* +Vi 

J m o Bessel function of the first kind of the m order 


) 


One can easily satisfy oneself that the expire ssion for q> 
remains unchanged in the case of opposite direction of rotation. 
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For the magnitude of the sound_ pros sure amplitude there results 

-i<= 

-A(R)e cos "5 






p E o 

55 p] 

s i2. 



3t 

2«r 

l 

i 

‘0 


+ B(R)e" ir, m 


J mn (3siR sin d) dR 


If one inserts 


A(R) dR = ~ 
n 


B(R) 6R = 


dM 

nR 


V 


xdo>l 


(oj^ «= circular frequency of the fundamental tone), P is 


P = 


mai^ 

2iter 


I R 0 

C e * €m cos d + *^~r ®e j mn (KR sin d) dR 
0 V dR cd x r 2 dR ; ® ' 


(9) 


The sound power may "be calculated according to the formula 


W 


f) » 2 

*^- 2 «r 2 sin d dd 

I r> 2P ° 

(Jo 


( 10 ) 


In order to execute in equations (8) and (9), respectively, the 
integration with respect to R, one would have to "know the distribution 
of the thrust end drag forces along the "blade \ it can "be determined far 
instance from aorodyinmic model tests. However, it is still possible 
to approximately evaluate the integrals if the aerodynamic data are 



12 


NACA TM No. 1195 


less detailed. Since 6 m and are small and the Bessel functions 

occurring in the integrand may "be regarded for the first harmonics as 
monotonically increasing with R it follows on the basis of the mean 
value theorem that: 


P = 


_ na>l 

2itcr 


-3 


•P coe d JjanflcRi sin d) + -^^J Bm (l£EU sin d) 

o>iR 2 2 


i 

- i 


( 11 ) 


Si and Eg are certain mean values. 

If the number of blades is small, one can, for the first 
harmonics, put Ri and Rq about equal to the radius E 0 of the 

circumference on which runs the point of application of the resultant 
thrust force of a single blade . In most cases it equals about 
0.7 - 0.75 Ro- 


KDn 

2«cr 


[< 


•P COS d + 


noM 

OJ3R 2 


Jmn(hE sin d) (R £ E 0 ) 


( 12 ) 


This contention shall be discussed for the oase of a two -blade 
propeller. 

In the integral 


1 Ro 
I/O 


~e~ iei J 2 (kR sin d) dR 


case may equate 


J 2 (KR sin d) a ~k e R 2 sin 2 d 

Q 
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One then obtains 


k 2 Bln 2 i 

“ 


)Bq 


iRO 


S e -Ui B 2 a . r a£ 


cffi 


8 


l 0 


dR 


* P 


k%i 2 sin^S 


8 


2 PJ 2 (1% sin 3) 


A comparison of the relations 


and 

r i Ka ^ - p%s 

shows that 




The upper limit is reached when the forces are distributed with 
the maximum of nonunif otmi ty , namely, concentrated at the ends of 
the blade. For the distributions occurring in practice, there is only- 
little difference between and S 0 » 
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In the Integral 



dM 

R2 dE 6 


iT,1 J 2 (lcR ein -3) dR 


the faotor 


J 2 (3sR ein $ ) 

_ 


is in first approximation fully independent of R and one can, 
therefore, substitute R <= Rjj then 



dM 

R2 dR® 




sin d) dR 


J 2 (1cRi sin d) 
R X 2 



dli -iTji 
dR 


dR 



J 2 (l® 2 _ sin d) 


Thus one obtains for the fundamental tone of a two -blade 
propeller the following expression for the Bound -pressure amplitude 


P s -SSL 0 

Sxor \ 


■J? cos d + 


sins) 


(13) 


(with R being slightly larger than R 0 ). 

For the second harmonio generally R^ 4 Rgj however, approximately 
one may assume also in this case 


P « 




C 


■P COS d + 


2CM 


sin d) 


2rtcr cojR^ ^4^ 

(ixi this case, the difference between R and R 0 is larger.) 


(HO 
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5. COMPARISON WHS THE EXPERIMENTAL RESULTS 


The theoretical results are now to he compared with the experimental 
data of Paris and Neap, for this purpose, the values P and M for 
the propeller used hy these investigators must he known; unfortunately, 
they are not given directly; however, they may he calculated approxi- 
mately bn the basis of other given data. Both investigators experimented 
with a two -blade propeller of diameter of 4.5 meters and a number of 
revolutions of 13*9 per second. The propeller operated under static 
conditions. If one assumes that the motor developed under static, 
conditions its full power - (600 hp), one obtains acoording to the 
formula (1) 


P » 1690 kg 

(under the assumption of i) a 0.75, a quite probable value). 
For the torque one obtains 

M « 515 kgm 


Using these numerical values, two directional characteristics 
were calculated for the fundamental acoording to the formula ( 13 ), 
corresponding to the two values R a 0 , 7 Rq and R a 0 . 75 Rq 
( cf. figs. 9 and 10 ). 

For the second harmonic a directional characteristic was 
calculated according to the formula (14) for the value S a 0 . 75 Rq 

(cf. fig. 11 ). 

A comparison with the experimental directional characteristics 
(cf . figs. 3 and 4) shows that the main characteristics of the 
experimental curves for the fundamental tone are well borne out by 
the theory. The agreement with Paris' results is particularly good. 

For the second harmonic the agreement with Kemp's curve (fig. 5) is 
less satisfactory. 

For the sound power results, acoording to the formula (10) by 
means of graphical and numerical integration, the values 34 W and 30 W 
for the fundamental corresponding to the values R = O.TRo and R * 0.75 Rq, 
and 6.5 W(R = 0.8Rq) and 4.7 W(B * 0.75 Kq) for the second harmonic. 
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The agreement with Kemp'e experimental values (17-8 W and 7*^ W) 
may he regarded as satisfactory if. one takes into consideration, on 
the one hand, the approximate character of the present calculation 
of P - and M, on the 9ther, the possible experimental errors, (in 
this connection, the noticeable quantitative discrepancy between 
Kemp's and Paris' curves for the fundamental should be pointed out.) 

Let it be noted that the formula (6) which forms the basis for 
the calculation of the sound field was derived from the equations of 
motion for small motions. This circumstance may lead to an under- 
estimation of the higher harmonics. Perhaps this is the reason for 
the increasing discrepancy between the theory and Kemp's experimental 
values for the higher harmonics. However, a final decision of this 
problem requires further and more accurate experimental data. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 1.- Calculated directional characteristic for the fundamental, 
tone of a two -blade propeller according to Lynam and Webb 
(hypothesis 1). 


f*90° 



Figure 2.- Calculated directional characteristic for the fundamental 
tone of a two -blade propeller according to Lynam and Webb 
(hypothesis 2). 


& 90 * 



Figure 3. - Experimental directional characteristic for the fundamental 

tone according to Paris. 
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Figure 4.- 


Figure 5.- 


#*9Q° 



Experimental directional characteristic for the fundamental 
tone according to Kemp. 

£* 90 ° 



Experimental directional characteristic for the second 
harmonic according to Kemp. 
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Figure 6.- Calculated directional characteristic for the fundamental 
tone according to Paris (combination of the two hypotheses by 
Lynam and Webb) . 



Figure 7. 


Figure 8. 
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Flight direction 

Figure 11.- Calculated directional characteristic for the second, 

harmonic. 



